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Abstract 

A method is proposed to create materials with a desired refraction coef- 
ficient, possibly negative one. The method consists of embedding into a 
given material small particles. Given no(x), the refraction coefficient of the 
original material in a bounded domain Del 3 , and a desired refraction 
coefficient n(x), one calculates the number N(x) of small particles, to be 
embedded in D around a point x £ D per unit volume of D, in order that 
the resulting new material has refraction coefficient n(x). 
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1 Introduction 



There is a growing interest to materials with the desired properties, in particular, 
with negative refraction coefficient (see [1] and references therein). In [2] the role 
of spatial dispersions is emphasized in explaining unusual properties of materials. 
In [3] the role of dispersion for wave propagation in solids is described. In jl] 
boundary-value problems in domains with complicated boundaries were studied. 
In [6], [7] wave scattering by small bodies of arbitrary shapes is studied and 
formulas for the S-matrix are obtained. In [5] a general method for creating 
materials with wave-focusing properties is proposed and justified. Our aim in 
this paper is to use a similar approach for creation of the materials with a desired 
refraction coefficient by embedding small particles into a given material with 
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known refraction coefficient no(x). The acoustic wave scattering by the given 
material is described by the Helmholtz equation 



[V 2 + k 2 n (x)]u = in M 3 , n (x) 



1 in D' := R 3 \ D, 

n (x) in D. 



Here k > is the wavenumber in D'. Equation (TjQ) can be written as the 
Schrodinger equation 

Lqu := [V 2 + k 2 — qo{x)]u = in R 3 , q := k 2 - k 2 n (x). (2) 

We assume k > fixed and do not show k- variable in go- Clearly, go = in D'. 
The scattering solution to ([2]) is uniquely defined by the radiation condition: 

Mo = e ifca - + A (/?,«) — + o(-), r := |x| ^ oo, /? := -. (3) 

Here a 6 5 2 is a given unit vector: the direction of the incident plane wave, S 2 
is the unit sphere in M 3 , A (j3,a) is the scattering amplitude, and (3 is the unit 
vector in the direction of the scattered wave. 

Assume that M small particles D m , 1 < m < M, are embedded into D. Small- 
ness means n ka <C 1, where a = 0,5max m diamD m , and n = m&x x( zD \no(x)\. 
On the boundary S m of D m an impedance boundary condition is satisfied: 

= CraU(s), S G S m , 1 < 7X1 < M, 

where is the unit normal to S m pointing out of D m . We assume that the surface 
S m is Lipschitz, and the Lipschitz constant does not depend on m, 1 < m < M. 
The scattering problem can now be stated as follows: 

M 

L Q u = in R 3 \ [J D m , u N = ( m u on S m , 1 < m < M, (4) 

u(x) = uq(x) + Am(P, a) r = Ixl — > oo, (3=—. (5) 

We prove that the solution to problem — (|5|) converges as M — > oo to the 
solution of the problem 

LU := [V 2 + k 2 - q(x)}U = in IR 3 , (6) 

U = e ika ' x + A{(3,a) ( -^- + o(^), r = \x\ -> oo, /3 = - , (7) 

where 

= q (x) +p(x), (8) 
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and give a formula for p(x). It turns out that p(x) can be made an arbitrary 
desired function by choosing the density of the number N(x) of the embedded 
particles around each point x G D and the impedances ( m properly Thus, q(x) 
can be made an arbitrary desired function. Therefore the refraction coefficient 

n(x) = 1 — k~ 2 q(x) = n (x) — k~ 2 p(x) (9) 

can be made arbitrary, in particular, negative. 

If no(x) is given and one wishes to create the material with the coefficient 
n(x), then one calculates 

p(x) = [n (x) — n(x)]k 2 , 

and embeds N(x) small particles per unit volume of D around each point x G D 
and chooses their impedances ( m so that the function p(x) is obtained for the new 
material. In Section [2] we give analytical formulas for N(x) and ( m and sufficient 
conditions for the convergence of the solution to (BJ - (J5j) to the solution of (J6]) - 
(IE|) as M — > oo in such a way that relations (|T3l)-(|T4l) hold. 

We also prove that the relative volume of the embedded particles is negligible. 
More precisely, if \D m \ the volume of D m , then 

Y M \D I 
lim ^ m=l 1 m| = 0. 

Af^oo D 

By \S m \ we denote the surface area of S m . We use an approximate formula 
for the electric capacitance of the perfect conductor with boundary S (see [6], 
p.26, formula (3.12)): 

(1"=^, J m := [ f (10) 



T / / k — t 

lib t/ o-m ,J Jm 1 

Note that 



C£ } = 0(a), \S m \ = 0(a 2 ), J m = 0(a 3 



By Cm the electric capacitance of a perfect conductor with the surface S m is 
denoted. We assume that 

n^ka <C 1, d^a, d :— mindist(-D m , Dj). (11) 

Let 

We assume throughout the paper that 

d = 0(a^), a = 0(i). 
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Let M — > oo and assume that the following limit exists: 

lim C^\Cm\S m \r l :=h{x). (13) 

M—>oo 
\x m -x\<d 

Here and below x m G D m is an arbitrary point in D m . Because D m is small, the 
choice of this point in D m is not important. Under the assumed relations between 
a and d one has lmiA/^oo § = 0. The limit (1T3T) exists if and only if ( m = 0(a _1 ), 
because \S m \ = 0(a 2 ) and Cm = 0(a). 

Denote by N m (x) the number of small particles per unit volume around a 
point x G D: Nn,i(x)dx = J2D m c£> 1 f° r an y subdomain D C D. 

The number of particles per unit volume is O(^) = (•=),. Therefore their 
relative volume is 0(^r) = 0(a 2 ) — > as M — > oo. On the other hand, the 
quantity NM(x)C m Q m , which has physical meaning of the average quantity C m ^ m 
per unit volume of D around point x, has a limit: 

lim N M (x)C mCm = -^L , lim N M (x)C^ := C(x). (14) 

M-+OC 1 + Ail X) M^oo 

\x m —x\<d \x m — x\<d 



The existence of the finite second limit in (1T4"1) is clear because Nm(x) = O(-) 

and Cf} = 0(a), and the existence of the first limit in (1T4"1) follows from formula 
(Tl3l) and from the second formula (HID. Our basic result is the formula: 



C(x) 

[wo(aO ~n(x)]fc 2 :=p(x) = , (15) 

1 + n[x) 

where C(x) is defined in ( TT41) and /i(x) is defined in ( fl3l) . 

Example 1. Suppose ( m = oo, so M\s m = 0' which corresponds to acoustically 
soft particles. Then h(x) = 0, p(x) = C(x). Assume that the small particles are 
balls of radius a. Then C^ = a, N M (x) = ^ , M = OQ). Since iV M (a;) > 
and Ct 0) > 0, then p(x) > 0, so one can create in this case only non-negative 
functions p(x). For any positive function p(x) one should embed N(x) = 
small acoustically soft balls of radius a per unit volume of D around each point 
x G D, and the resulting material will have n(x) = rio(x) — k~ 2 p(x). In particular, 
n(x) < if p(x) > k 2 rio(x). 

Example 2. Choose an arbitrary function p(x) = p\(x) + ip2(x), p2(x) < 0. 
The condition pi < guarantees uniqueness of the solution to problem (jHJ)- 
(JTj) with q(x) = qo(x) +p(x). Physically this condition means that the medium, 
corresponding to n(x) = l—k~ 2 q(x) has nonnegative absorption. Let the particles 
be balls of radius a and ( m = ( m ( x) — 4 L nah(x) ' where h(x) is an arbitrary function 
at the moment. This function is fixed later. Then formula (flBl holds because 
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\S m \ — Ana 2 . Choose N = N(x) and h(x) — h\ + ih 2 from the first equation 
using (H2J: 



iVa iVa(l + h x - ih 2 ) 

P1+P2- 



Thus, 



l + h{x) {l + hx) 2 + hl ' 
Na{l + hx) Nah 2 

We have three functions: N = N(x) > 0, hx and h 2 , to satisfy two equations 
(JTBT) . This can be done by infinitely many ways. For instance, one can take 
hx — 0, hi — — ^ and iV = a -1 pi(l + ^f). Thus, to get the material with 
the desired n(x) = n (x) — k~ 2 p(x), where p(x) = px(x) + ip 2 (x), one embeds 
N(x) = a~ 1 {p\ +p 2 )/px small balls of radius a per unit volume around each 
point x and chooses the impedance ( m (x) = (47ra/i(x)) , where h = hx + ih 2 , 
h 2 = -p 2 /px, hx = 0. 

2 Derivation of the results. 

We seek the unique solution to (j4j) - ((Sj) of the form 



M 

u = u + 

m=l 



J2 f G(x,t)a m (t)dt (17) 

m=l J Sm 
M M r 

uo + y^ j G(x,x m )Q m + y~' j / [G(x,t)-G(p 

m=l m=l J 



- G(x,x m )\a m dt. 



Here L0G1 = 5(x — y) in IR 3 , G satisfies the radiation condition, a m are to be 
chosen so that the boundary condition (TjJ is satisfied, Q m := f s o m dt, x m G D m . 
In the generic case Q m 7^ one can neglect the last term in ffTTl) compared 
with the preceding term if \x — x m \ > d ^> a for all m. Indeed, under this 
assumption one has \G(x,t) — G(x,x m )\ < \V y G(x,y) ■ (t — x m )\ = O(^) 1, 
where y := x m + r(t — x m ), < r < 1, is a 'middle point'. Thus, the third term 
on the right side of (|T7|) is \Q m \) \Q m \, where we also assume that \Q m \ = 
0(J S \cr m \dt). We will see that this assumption is justified. For example, if 
u\g =0, then a m does not change sign on S m . Thus, generically one can write 



M 

il 

m=l 



u a( x ) + G (x,x m )Q m , \x - x m \ > d > a, (18) 



with the error 0(%), The choice of x m G D m does not matter because a is small. 
One may assume that D m are convex and take x m at the gravity center of D m . 
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The functions G(x,y) and uq(x) are known because uq{x) is known. Let us 
derive an equation for finding Q m . If Q m are found then the scattering problem 
(HJ) - (EJ) is solved by formula (JIB!) for any x away from an immediate neighborhood 
of the small particles. To derive an equation for Q m we need some preparations. 
The function G(x, y) solves the equation: 

r e %k\x-y\ 

G(x,y) = g(x,y) - g(x, z)q(z) G(z,y)dz } g(x,y) := — (19) 

Jd An\x-y\ 

One can easily prove that 

G{x, y) = g{x, y) [l +0{\x- y\)] = g {x, y) [l + 0{\x- y\)] , \x - y\ -+ 0, (20) 
where go(x, y) = (47r|x — Let 

dg (s,t) 



aJt)dt. (21) 



T 3 a 3 := ! G(a,t) affldt, Vi = 2 / <)X 

J Sj ^ 

It is known (|6j, p. 91 ) that 

^* = - /„ jW *, gffi2) = ^-^ , (22) 

where Aj(k) is the operator similar to (1211) with g(s, £) in place of go(s, t), N s := iV 
is the outer normal to Sj at the point s G Sj. On the surface Sj we have 

M 

u = u e (s) +Tj(Tj, u e := u + 2j G(s,x m )Q m . (23) 

Using boundary condition (J4J) and formulas (1221 . (|23|) . one gets 

Mejv (s) - 0« e (s) + Aj<7j ~ aj - QTjCTj = 0. (24) 
Integrate (|24|) over S^, use (1221) and get: 

Qj f = / u eN (s)ds — (j / u e (s)ds — Q / TjVjds. (25) 

J S'i J On J O n 



One has 



/ u eN ds = I Au e dx = 0{k 2 a?) 1 I u e ds = u e {x 3 )\Sj\ = 0{ 

J Si JDi JS-i 



a 2 ), 



where the smallness of Dj and the fact that u e and its two derivatives are bounded 
on Sj were used. Since fca C 1 we can neglect the first integral in ( |25l) compared 
with the second. Furthermore 



1 ■= T j°j ds = I dtaj{t) j 

J S'i J S -i *-* 5* -j 



ds 



4tt\s - 1\ 
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We replace the last integral by its mean value 



dt 



ds Jj 



Jsj JSj 47r|s — t\ 47r|Sj| 



Thus, I = , and ( 1251) yields: 



* = ~ ■ iMm^ u - (Xi) ■ (26) 



We have replaced u e (s) by u e (xj) because \xj — s\ < 2a and a is small while u e (x) 
is continuous in a neighborhood of Xj. Using (flQj) we rewrite ff26|) as 

Qi = ~Cf[l + Cf {QlSjir 1 }- 1 ^) := -C Kj u e (xj). (27) 
Thus, (JTHJ) can be written as: 

M 

u(x) = u (x) - G(x, x m )C mCim u(x m ), \x - x m \ > d^ a. (28) 

m=l 

We have replaced -u e (x m ) by u(x m ) under the sign of the sum in ( 1281) because at 
the points x which are away from small particles one has u e {x) = u(x) with the 
error 0(2). Formulas flU - ([HD allow one to pass to the limit M — > oo in (1281) 
and get 

W(x)=uo(ar)- / G(x,y) p(y)U(y)dy, (29) 
in 

where p(x) is defined in (1151) . Applying to (1291) the operator Lq, defined in (J2J), 
and using the relation L G(x, y) = — S(x — y) yields equation (jSJ) with g defined 
in ([S]). The radiation condition for W is satisfied: 

A(/3, a) = Ao(/3, a) + Ai(/3, a), (30) 

where 

A x (/3,a) = lim A M (P,a) = ~ [ u (y,-j3)p(y)U(y)dy. (31) 

M->oo 47T Jj) 

Here we have used a result from [5]: 

= ——u {y, -/?) + o(— r), |x| -»• oo, -p-r = j3. (32) 
47r|x| N |x| 

In our derivations it was assumed that ( m 7^ 0. If ( m = for all m, that is, the 
small particles are acoustically hard, then Q m = in the first order with respect 
to ka. One can show that in this case Q m = 0(k 2 a 3 ), and that the last sum 
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in (ITTj) is of the same order of magnitude as the preceding sum. Consequently, 
the theory in this case is quite different: the effective field in the medium is not 
described by equation (!29j) . which is equivalent to a local equation (jHJ). In fact, 
the effective field in this case is described by an integrodifferential equation which 
is not equivalent to a local differential equation. 

Let us explain the relation Q m = 0(k 2 a 3 ), mentioned above. Write (T24T) with 
Q — 0, integrate over Sj and use the first formula ff22l) to get 

Q j — / u eN ds = / Au e dx = 0(k 2 a 3 ). 
Js j Jd 3 
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